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Consequences of strong fluctuations on high-energy QCD evolution 

C. MarquetQ R. PeschanskiQ and G. Soyezfl 
Service de physique theorique, CEA/Saclay , 91191 Gif-sur-Yvette cedex, France' 

We investigate the behaviour of the QCD evolution towards high-energy, in the diffusive approx¬ 
imation, in the limit where the fluctuation contribution is large. Our solution for the equivalent 
stochastic Fisher equation predicts the amplitude as well as the whole set of correlators in the 
strong noise limit. The speed of the front and the diffusion coefficient are obtained. We analyse the 
consequences on high-energy evolution in QCD. 


1. The quest for the perturbative high-energy limit of 
QCD has been the subject of many efforts. It is now well 
accepted that, due to the strong rise of the amplitude 
predicted by the linear Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation Q, one has to include saturation effects 
in order to describe high parton densities and recover uni- 
tarity. In the large- N c limit and in the mean-field approx¬ 
imation, we are led to consider the Balitsky-Kovchegov 
(BK) equation J2]. It has the nice property 0 to be 
mapped, in the diffusive approximation, onto the Fisher- 
Kolmogorov-Petrovsky-Piscounov (F-KPP) equation jj| 
which has been widely studied in statistical physics and is 
known to admit traveling waves as asymptotic solutions, 
translating into geometric scaling in QCD @. 

It has recently been proven jj| Q that fluctuation ef¬ 
fects have important consequences on the approach to 
saturation. Practically, the resulting evolution equation, 
after a coarse-graining approximation Q, takes the form 
of a Langevin equation. It is formally equivalent to the 
BK equation supplemented with a noise term 

dyT(L,Y) = a X (-d L )T(L,Y) - aT 2 (L,Y) 

+ a^KalT(L,Y)v(L,Y) (1) 

where T is the event-by-event scattering amplitude, Y 
is the rapidity, L = log(fc 2 /fcg) with k the transverse 
momentum and kg some arbitrary reference scale. x( 7 ) = 
2^(1) — ipi'y) — ip( 1 — 7) is the BFKL kernel, n is a fudge 
factor and the noise is(L, Y) satisfies {v) = 0 and 

(v(L, Y)u(L', Y')) = —5(aY - aY')5(L - L'). (2) 

7r 

To obtain equation 0, we have worked with impact- 
parameter-independent amplitudes, for which the orig¬ 
inal non-local, off-diagonal noise term takes HI the 
form given by equation 0 . Physically, in addition to the 
BK saturation effects coming from pomeron merging in 
the target, one also takes into account pomeron splitting 
in the target. Hence, by combinations of splittings and 


mergings, pomeron loops are formed. The extra factor 
a 2 in the fluctuation term indicates that it is dominant 
when T ~ a 2 i.e. in the dilute limit. 

In the same line that the BK equation is equivalent to 
the F-KPP equation in the diffusive approximation, the 
Langevin problem corresponds to the stochastic F-KPP 
(sFKPP) equation 0. To be more precise, let us expand 
the BFKL kernel to second order around 70 

x(7) = Xo+Xo(7-7o) + ^Xo(7-7o) 2 = A 0 +A lX +A 2X 2 . 

(3) 

This approximation has proven its ability to exhibit the 
main properties of the solutions of equation 0 in the 
limit na 2 <C 1. Unless specified, we shall keep this ap¬ 
proximation throughout this paper, leaving the general 
case for further studies. 

If we introduce time and space variables as follows 

T 

t = aY, x = L — A\aY and u{x,t) = 

Ao 

equation 0 gets mapped onto the sFKPP equation 1 
dtu = Dd 2 u + Au(l — u) + E\J u{\ — u)rj(x , t), (4) 


with 


D = A 2 , A = Aq, e 2 
(ri(x,t)r)(x' ,t')) = S(t- 


2 na 2 

= Tv 

t')5(x — x'). 


At present stage, most of the analytical analysis were 
performed in the limit where the fluctuations are asymp¬ 
totically small (the correction being logarithmic |1(1 . 
it may require a strong coupling constant as small as 
07 1 $ 1CD 10 ), in which case the relevant quantities, e.g. 
the speed of the wave, can be expanded around the F- 
KPP solution. In this analysis, the main effects of the 
noise in the sFKPP equation are, first, to lead to a de¬ 
crease of the speed of the traveling front and, second, to 
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We have introduced and extra factor \/l — u in the noise term. 
The effects of the noise being important in the dilute tail, this 
modification is not expected to change physical results. In addi¬ 
tion, in q. the fluctuation contribution is only under control in 
the dilute regime. 
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produce a diffusion of the fronts for each realisation of 
the noise resulting in violations of geometric scaling for 
the average amplitude. 

The numerical studies performed so far show that these 
effects (decrease of the speed and diffusion of the events) 
are amplified when the fluctuation coefficient becomes 
more important. There has been large efforts made to 
improve the analytical understanding for arbitrary values 
of the noise stren gth but many approaches appear to be 
model-dependent [ill . 

In this paper, we consider the problem of fluctuations 
through a complementary approach, namely the limit of 
a strong noise. This limit is tractable with the help 
of a duality property of the sFKPP equation |§j. The 
strong-noise limit then gets related to a coalescence pro¬ 
cess which can be solved exactly El- 

Using these tools from statistical physics, we are able 
to compute the event-averaged amplitude as well as the 
higher-order correlators and obtain predictions for the 
speed of the wavefront as well as for the diffusion coeffi¬ 
cient in the limit of strong fluctuations. This knowledge 
of the strong-noise limit, together with the weak-noise 
results, can help in further analytical understanding of 
the QCD evolution in the presence of fluctuations. 

2. Let us now summarise the tools from statistical 
physics we need for our studies. Our starting point is the 
duality relation between the sFKPP Langevin equation 
and the reaction-diffusion process. This duality will allow 
us to relate the strong noise problem to the coalescence 
problem which we shall solve. 

We consider on the one hand amplitudes evolving ac¬ 
cording to the sFKPP equation © and, on the other 
hand, the reaction-diffusion process of a one-dimensional 
population A on a lattice of spacing h: at each site, one 
can have particle creation or recombination, and particles 
can diffuse to a neighbouring site 


Therefore, to obtain information on the evolution of 
the average amplitudes for the sFKPP equation, we shall 
study the dual particle system and then use relation ©• 
The limit we are interested in is the strong noise limit 
(large e). This corresponds to heavy saturation in the 
particle system (projectile), i.e. two particles at the same 
lattice site automatically recombine into a single one. 
This limit is often referred to as the diffusion-controlled 
limit. When A is rescaled to give a constant, small, ratio 
A/e 2 , we can alternatively study the coalescence model. 
In this model, one can have at most one particle per lat¬ 
tice site. One particle can diffuse to the neighbouring site 
at rate D/h 2 or give birth to a new one in a neighbouring 
site at rate ui/h (with uj = 2D A/e 2 as we shall see later). 

This system has been studied 0 and is fully solvable 
using the method of inter-particle probability distribution 
function. The main idea is to introduce E(x,y;t ) as the 
probability that sites between x and y > x included are 
empty at time t. One obtains that, due to diffusion and 
creation, E satisfies the following differential equation 

d t E = {D (,9 2 + d 2 ) + u; (d y - d x )} E. (7) 

with the boundary condition limj,-*^ E(x , y\ t) = 1. 

The particle density can be obtained from the deriva¬ 
tive of E: 

N(x,t) = d y E(x,y;t)\ y ^ x . 

Remarkably enough, the evolution equation 0 is linear. 
It can be solved exactly [3 and, for a given initial con¬ 
dition E(x,y,t), introducing the dimensionless variables 

u> , . u> . . 8 u> 2 

£, = jj(x + y), (=p{y-x) and r = — t, 

one finds 


Ai^Ai + Ai, A+A^A, and A°-^ A l±1 

( 5 ) 

where A, designs a particle at lattice site i. 

One shows Q that the particle system and the ampli¬ 
tude in the sFKPP equation are related through a dual¬ 
ity relation which, in the continuum limit h —> 0, can be 
written 

/n [i - «(*, or 0 - 01 )=(n n - «<*. . 

( 6 ) 

where N(x, t ) is the particle density in the reaction- 
diffusion system. Physically, this duality equation means 
that, if one wants to obtain the scattering amplitude at 
rapidity t = aY , one can either keep the target fixed and 
evolve the projectile wavefunction considered as a parti¬ 
cle system, or fix the projectile and consider evolution of 
scattering amplitudes off the target. One knows 013 
that splitting in the projectile leads to linear growth and 
saturation in the target while merging in the projectile 
corresponds to fluctuations in the target. 


E{x,y;t) 








( 8 ) 


where the Green function G(£, (, r) is given by 


_L e -(£-0 2 U- e -(C-C')/2 

7TT 


e -(C-C') 2 /r _ e -(C-K') 2 /r 


Before considering the solution of this system in the 
context of the duality relation, let us give the relation be¬ 
tween the parameters A and e 2 of the initial system with 
ui and D in the coalescence model. The trick is to require 
that both systems have the same equilibrium density. For 
the coalescence model, one notice that exp(— fj{y — x)) is 
a time-independent solution leading to a particle density 
N eq = tu/D. In the case of process ©, at equilibrium, 
diffusion does not play any role and we have to find equi¬ 
librium at each site between creation end annihilation. 
This is achieved when N eq = 2 A/e 2 . Hence, one has 
u> = 2 D\/e 2 . 
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3. Let us now put together the results from duality 
and coalescence and derive the sFKPP solution. 

By carefully choosing the initial condition, equation 
© simplifies. If one starts with one particle at po¬ 
sition x i.e. N(x, 0) = 8(x — x), then the l.h.s. of 
© becomes simply 1 — (u(x,t)). By starting with k 
particles at position x\,...,x k , one similarly obtains 
([1 - u(xi,t )\... [1 - u(x k ,t)]). 

In addition, let us start with a step function for the am¬ 
plitude zt(x, 0) = 6{—x). Then ^TIj; [1 — u i x > 0)] JV ^ :,t ^ 
is the probability that, in the particle process, all sites 
x < 0 are empty. For the case of the strong noise i.e. 


when the particle system is the coalescence model, this 
probability is directly obtained in terms of the density E 
and the duality relation becomes 

( u(x,t)) = 1 - E x (-oo,0;t), 
with the initial condition 

E x (x,y\ 0) = 1 - 6(x - x)6(y - x). 

Inserting this initial condition inside the general solu¬ 
tion ©, we get after a bit of algebra 


J 


E^y-t) = I|erfc(^-^^)+erfc 


x — y—2u>t 


— erfc 


x — y—ujt 


1 


+ -e d 




2 — erfc 


y—x—2ujt\ 

V8Dt ) 


V8Dt J \ 2 y/Dt 


1 —-erfc 


— x—ujt 

2VDt 


x — X + Ult\ 

2y/Dt ) _ 


1 —j-erfc f x ~y~^ Lot 

2 l 2 Vm 


where erfc(x) is the complementary error function. The We need to use this relation to obtain the correlations of 
limit x —-> — oo, y —> 0 in this expression gives u instead of 1 — u. This is obtained as follows 


<«<*,(» = Ierfc( f^) (9) 

(i-ul ) 2 

= — , / dxdix— x)e 4D * . 

2sfD^tJ- 00 

This results calls for comments. First, it corresponds 
to a wave traveling at an average speed 


2DX 



This result confirms the decrease of the velocity but 
contrasts with the speed obtained in the weak noise 
limit by perturbative analysis around the F-KPP speed 
~ 2VDX - n 2 VDX\ log(e)|” 2 . 

The expression © shows that the amplitude could be 
obtained from a superposition of step f uncti ons around 
x = u>t with a Gaussian noise of width \j2Dt. The inter¬ 
esting point here lies in the dispersion coefficient: in the 
weak-noise analysis, it behaves like |log(e:)|~ 3 . We have 
proven that this rise goes to 2D when the noise becomes 
strong. 

In addition, one can probe the correlators of the ampli¬ 
tude by starting with an initial condition with particles 
at positions x m - ln = X\ < ■ ■ • < x k = £ max : 

k 

E(x, y\ 0) = 1 - ^2 - Xi-i)0(xi - x)6(y - Xi), 

i=i 

with, formally, Xq = — oo. Following the same lines as 
above, one gets 

([1 - u(xi,t)} ... [1 - u(x k ,t)]) = 1 — 7 ^erfc ) • 


{u(xi,t). ..u(x k ,t)) 

= ({!-[!- u(xi , t)]} ... {1 - [1 - u(x k , t)]}) 

AC{l,...,fc} \ieA 

= E ( 1 - U ( X min(A) , t)) 

AC{l,...,Jfc} 

where jj^4 is the cardinal of the set A and min(H) is its 
minimum. The sum can be reordered to give 

(u(x!,t) ■ ..u(x k ,t)) 

k 

1=1 AC{j+l,...,fe} 

j =1 n—0 ' 

In the last expression, only the term with j = k survives 
hence (with cc max = x k ) 


(u(xi,t).. .u(x k ,t)) = -erfc 


x k — U)t 

2Vm 


(u(x k ,t)) ■ 


This simple result can again be seen as a superpo¬ 
sition of step functions with a Gaussian dispersion. 
u(xi).. .u{x k ) is nonzero provided u does not vanishes 
at the position with largest coordinate (x k ). Hence, the 
whole dynamics is the same as if only one particle were 
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FIG. 1: Comparison between the F-KPP and sFKPP wave- 
fronts at t = 0, 5,10,15 and 20. Solid curve: numerical simu¬ 
lation of the F-KPP equation. Dashed curve: analytic result 
for e 2 = 1.5. Dotted curve: analytic result for e 2 = 3. 

lying at position x k in the initial condition 2 . 

In order to illustrate the effect of the noise, we have 
plotted in figure ^ the time evolution of the wavefront 
(u(x, t)) for numerical simulations of the F-KPP equation 
and for our solution m (e 2 = 1.5 and 3). In each case, 
the initial condition is a step function. We clearly see that 
the effect of the strong noise is to slow down considerably 
the front and add significant distortion. 

4. Coming back to QCD variables, we find (assuming 
without loss of generality L\ < ■ ■ ■ < L n ) 

(T(L 1 ,Y)...T(L k ,Y)) = A *- 1 (T(L k ,Y)) (10) 

= 4 ertc ( U - caY \ 

2 VV2 D^SYj 

with the speed of the traveling front and the diffusion 
coefficient 

n T A A 2 

c = A\ -|-and = 2 A 2 . (11) 

naj 

Let us discuss the physical interpretation of these re¬ 
sults. We start by equation El. It is remarkable that 
the strong-fluctuation limit gives, as an analytic solution, 
the asymptotic result inferred in previous studies nm 
It proves the universal feature that high-energy scattering 
amplitudes are described by a superposition of Heavyside 
functions with Gaussian dispersion. We confirm analyti¬ 
cally that the correlators are driven by the amplitude of 
the largest momentum in the process. This is the main 
result of this paper. 


2 Of course, the same argument holds for the product of 1 — u(xt ) 
and the fronts around x\. 


Equation CD relates the parameters of the amplitude 
e.g. the average speed c and the diffusion constant Ddiff 
to the parameters Ao, A\ and A 2 obtained from the ex¬ 
pansion 0 of the BFKL kernel. If one performs this ex¬ 
pansion choosing 70 =\ or 70 = 7 c ~ 0.6275, as used in 
the weak-noise limit, the value of A\ turns out to be neg¬ 
ative. This would lead to a negative speed which seems 
unphysical in QCD. The way out of this inconsistency is 
to insert the solution m directly into the exact evolu¬ 
tion for (T(L, Y)}. It has been shown El that it results, 
as expected, into a vanishing speed. The determination 
of c through the evolution equation would depend on the 
corrections to the error function, which disappear in the 
strong-noise limit. 

Let us sketch a heuristic argument indicating that a 
physically meaningful behaviour of these parameters can 
be obtained in the limit 70 small. Indeed, this is sug¬ 
gested by the fact that, in the strong-noise limit, each 
front is approaching a Heavyside function © which sug¬ 
gests to perform the kernel expansion near 70 = 0. Con¬ 
sidering 70 small, one has Aq = SJq 1 , Ai = — 37,)“ 2 and 
A 2 = Jq , leading to 

7o \7o™s / 

If one requires c —> 0, one has to choose 



where o(x) denotes a function falling to zero faster than 
x. The physical parameters, in the strong-noise limit, are 
then 

O kcA 

c - 0 and £> diff = —X ( 12 ) 

37T 

It is interesting to notice that, if we choose 70 to satisfy 
the only requirement that c —> 0 , the diffusion coeffi¬ 
cient is entirely determined, independently of the way 
this limit is achieved. 

On a more general ground, it is worth noting that the 
strong-noise development should not be considered as a 
strong coupling expansion since the initial equation 0 
is derived in a perturbative framework. Indeed, the gen¬ 
uine expansion parameter in the strong noise limit is k « 2 
where n could reach high values in the physical domain of 
interest fl6§ . Thus, even in the perturbative regime, the 
noise parameter na 2 s may be large. These predictions, 
together with the weak noise results, enclose the physi¬ 
cal domain of na 2 s . The knowledge of both these limits, 
together with information from numerical estimates [l7| , 
can lead to a better understanding of the physics of fluc¬ 
tuations. An expansion in 1 /{not 2 ) could also lead to 
faster convergence than the logarithmic weak-noise ex¬ 
pansion pH |. 

5. Let us summarise our results. We have used 
the duality relation between the amplitude given by the 
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stochastic FKPP equation and the particle densities in 
a reaction-diffusion process. This duality is physically 
similar to the projectile-target duality noticed recently 
in high-energy QCD when both saturation and fluctu¬ 
ation effects are taken into account. The saturation in 
the target is related to splitting in the projectile while 
fluctuations are mapped to recombination in the particle 
system. 

In the case of large fluctuations i.e. strong noise in 
the sFKPP equation, the corresponding particle system 
can be described as a coalescence problem. This process 
can be solved exactly using the inter-particle probability 
distribution function. We use this to compute the average 
value of the amplitude as well as the correlators. 

The main result of our analysis is the analytic deriva¬ 
tion of the average scattering amplitude as a universal 
error function m which also determines higher-order 
correlators. This picture proposed in previous studies is 
thus confirmed and shows that the results obtained in the 
limit of strong fluctutions possess a physical meaning. 

The fact that the correlators display the same be¬ 
haviour as the amplitude itself is physically interesting. 
This is obtained through a superposition of event-by¬ 


event amplitudes which are 0 or 1. The dominant contri¬ 
bution to the scattering process comes when all individ¬ 
ual scattering are 1. This picture of a black and white 
target gives rise to new scaling laws for Deep Inelastic 
Scattering EH- 

Since the physically acceptable values of the strong 
coupling seem to lie in between the strong and the weak 
noise limits, a knowledge of both approaches is useful. 

As an outlook, it would be interesting to extend our 
formalism beyond the diffusive approximation and/or to 
modified evolution kernels [l8|. This may allow a bet¬ 
ter determination of the speed and dispersion coefficient, 
given by m in the diffusive approximation. Also, a per¬ 
turbative approach starting from the strong noise limit 
could prove useful. These questions certainly deserve fur¬ 
ther studies. 
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